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Abstract. We consider the one-dimensional XY quantum spin chain in a transverse 
magnetic field. We are interested in the Renyi entropy of a block of L neighboring 
spins at zero temperature on an infinite lattice. The Renyi entropy is essentially the 
trace of some power a of the density matrix of the block. We calculate the asymptotic 
for L — > cxD analytically in terms of Klein's elliptic A - function. We study the limiting 
entropy as a function of its parameter a. We show that up to the trivial addition 
terms and multiplicative factors, and after a proper re-scaling, the Renyi entropy is 
an automorphic function with respect to a certain subgroup of the modular group; 
moreover, the subgroup depends on whether the magnetic field is above or below its 
critical value. Using this fact, we derive the transformation properties of the Renyi 
entropy under the map a —^ a^^ and show that the entropy becomes an elementary 
function of the magnetic field and the anisotropy when a is a integer power of 2, this 
includes the purity trp^. We also analyze the behavior of the entropy as a and 
CXI and at the critical magnetic field and in the isotropic limit [XX model]. 
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1. Introduction 

Entanglement is a resource for quantum control [1]. It is necessary for building quantum 
computers. Different measures of entanglement are used in the literature. For pure 
systems [considered here] the von Neumann entropy of a subsystem is the most popular 
measure [21 [Sj HI [5l [6l [7]. The subsystem is a large block of spins in the unique 
ground state of a spin Hamiltonian. In this paper we evaluate the Renyi entropy of the 
subsystem. The Renyi entropy was discovered in information theory P [HI [101 [HI [12] , it 
is essentially the trace of a power of the density matrix. For physics the Renyi entropy 
is important, because once we know the value of the trace of every power of the density 
matrix then we can reconstruct its whole spectrum. 

The physical system we consider is the anisotropic XY model in a transverse 
magnetic field and the entropy we are interested in is the one of a block of L neighboring 
spins at zero temperature and of an infinite system. The Hamiltonian for this model 
can be written as 

oo 

H = -Y,i^ + 7)«+i + (1 - + ha:, (1) 

^^=— oo 

Here < 7 is the anisotropy parameter; cr^, cr^ and cr^ are the Pauli matrices and 
> h is the magnetic field. The model was solved in [T3l [HJ [151 HH]- We are going to 
calculate the bipartite block entropy of the ground state \GS) of the system. 

The XY model can be mapped exactly into a system of free fermions with spectrum 
given by 

efc = A^{cosk-h/2f + sin^ k . (2) 

We can read the phase diagram of the model from its spectrum and identify that it 
is critical for 7 = 0, /i < 2 (corresponding to the isotropic XY model, or XX model) 
and at the critical magnetic field h = he = 2. For h = /i/(7) = 2a/1 — 7^ (factorizing 
field) the ground state can be written as a product state, as it was found in [T7j, and 
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Figure 1. Phase diagram of the anisotropic XY model in a constant magnetic field 
(only 7 > and h > shown). The three cases 2, lA, 1b, considered in this paper, 
are clearly marked. The critical phases (7 = 0, ft. < 2 and h = 2) are drawn in bold 
lines (red, online). The boundary between cases lA and 1b, where the ground state 
is given by two degenerate product states, is shown as a dotted line (blue, online). 
The Ising case (7 = 1) is also indicated, as a dashed line. 
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is doubly degenerate: 

\GS,)= n [cos(^^)| Tn) + sin(e)| U] , 
nglattice 

\GS2)= n [cos(^^)|Tn)-sin(^^)|U] , (3) 
nglattice 

where cos^(26') = (1 — 7)/(l + 7). Off this line, the ground state of the model \GS) is 
in continuity with the state 

\GS)h=hf{'Y) = \GSi) + \GS2) ■ (4) 

The line h = hf{'~f) is not a phase transition, but the entropy has a weak singularity 
across it, since its derivative, although finite, is discontinuous. In Fig. [T]we show the 
phase diagram of the XY model and mark the three regions where we calculate the 
different expressions of the entropy. 

We shall calculate the entropy of a block of L neighboring spins (a subsystem) of 
the ground state \GS) as a measure of the entanglement between this block and the rest 
of the chain. We treat the whole chain as a binary system \GS) = \A^B). We denote 
this block of L neighboring spins by subsystem A and the rest of the chain by subsystem 
B. The density matrix of the ground state can be denoted by pab = \GS){GS\. The 
reduced density matrix of subsystem A is pA = TrBipAs)- Then, the von Neumann 
entropy S{pa) and the Renyi entropy Sa{pA) of the block of spins can be evaluated by 
the expression 

S{pa) = - rr(pAlnpA), (5) 

SJpA = ^—lnTr{pVj, a 1 and a > 0. (6) 
1 — a 

Here the power a is a parameter. When evaluated for 1-dimensional critical theories, 
these entropies diverge logarithmically with the size of the block, while they saturate 
to a constant in the presence of a gap |18j . 

For the isotropic version of the XY model 7 = we evaluated Renyi entropy of 
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a large block of spins in [IS]. The von Neumann entropy of the block in the XY 
model was calculated in [201 [221 123]. The methods of Toeplitz determinants 
[2H [25l [261 [271 [281 [29] . as well as the techniques based on integrable Fredholm 
operators [30l [311 [22], have been used for the evaluation of the von Neumann entropy 
of this model [HI [33]. 

In this paper we evaluate the Renyi entropy, which is the natural generalization 
of the von Neumann entropy [8]. When the Renyi entropy turns into the von 

Neumann entropy. 

2. Renyi Entropy. 

The von Neumann Entropy of the block of spins has been calculated in [33] and [7]. 
We shall use the same notations and introduce an elliptic parameter: 



^{h/2f + 7^-1/7, 



Case la: 4(1 - 7^) < /i^ < 4; 



k=l ^(l-/i2/4-72)/(l-/i2/4) , Case lb: /i^ < 4(1 -72); (7) 



^ 7/ v/(V2P + 72-l, 



Case 2:h>2. 



We shall also use the complete elliptic integral of the first kind 
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and the modulus 



To = i{k')/m, 



k' = Vi^, 



(9) 



as well as: 



e = 



q = e 



^-nlik')/lik) 



(10) 



We will need the following identities as well |34j : 




(11) 
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Now let us start the evaluation of the Renyi entropy of a block of L neighboring 



spins. It can be represented \i9\ as 

L 



(13) 



where the numbers 



±ii/fc, k = l,...,L 



are the eigenvalues of a certain block Toeplitz matrix. In [33] it is shown that in the 
large L limit the eigenvalues z/2m and i'2m+i merge to the number defined below in 
eqm, 

Hence, the Renyi entropy in the large L limit can be identified with the convergent 
series, 



m — — rv~i I- ^ ' ^ 



with 



Am = tanh I m H — | ttto. 



(14) 



(15) 



The summation of the series can be done following the same approach as in the case 
of the von Neuman entropy (cf. ^). 



2.1. h>2 

With e = TTTo, we have 

1 + Xm = 2 



1 



I _|_ g-(l+2m)e 



-(l+2m)e 



1 -^m - 2^ _^^_^^_^2m)£' 



(16) 

(17) 



Renyi Entropy of the XY Spin Chain 



7 



Then, the entropy is: 



S 



1 — a ^-^ 



m=—oo 
oo 



1 - a ^ 



m=0 

oo 



X -)- g-(l+2m)e 
1 

X _)_ g-(H-2m)e 



+ 



-(l+2m)e 



+ 



X -|- g-(l+2m)e 
-(l+2m)e ^ 



X -|- g-(l+2m)e 



^Vln 

1-a ^„ 



m=0 
oo 



l + e" 



-C((l+2m)e 



:i + e-(i+2H0 



= ^— 5^ In [1 + e-"(^+2n^).] _ ^ ^ In [1 + e-^^+^m).] 

m=U m=0 

Summing the second term is straightforward, using ffTTj) : 

oo oo 



m=0 



m,=0 



1 a 
T21-a 



In g + In 



16 



/2 



where, as usual. 



g = g-7r/(fe')/^(fc)_ 



(19) 



(20) 



In order to sum up the first term we notice that identity ( fTTi) can be interpreted 
as the evaluation of the product in the left hand side in terms of the function k = k{q) 
defined implicitly by equation (1201) . A fundamental fact of the theory of elliptic 
functions is that the function k{q) admits an explicit representation in terms of the 
theta-constants. Indeed, the following formulae take place (see e.g. [M]): 



k{q) 



k\q) 



^1(0, g) 



21] 



^1(0, g)' 91(0, q)' 

where 9j{z\q), j = 1,2,3,4 are the Jacobi theta-functions. We remind (see again ^ 

that the theta functions are defined for any |g| < 1 by the following Fourier series 



n=—oo 



eoiz 



n=—oo 



(22) 
(23) 
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2izn 



n=—oo 
oo 



(24) 
(25) 



ra=— oo 

In particular, it follows that the functions 

k'{q) and q'^''^k{q) , (26) 

are analytic on the unit disc |g| < 1. It is also worth mentioning the classical formula 
for the integral /(fc), 

/(fc) = 1^32(0, g). (27) 

Put now, 

K := fc(g") , (28) 

where q is the q - parameter corresponding via (1201) to the original elliptic parameter 
k from (171). Then for the first term in ffTSl) we will have 



9 _ _ 9 °° 

^ In [1 + e-"(^-^-)^] = ^ In n (1 + (^")^'""^) 

m=0 m=0 



1 1 

121-a 



a In g + In 



16 
klk'^ 



(29) 



Substituting this expression together with ([T9|) into ([T8|) . we arrive at the equation. 



5, 



1 1 , / 16 
In 



12 1 - a V klk';^ ) 121 -a 



1 a , / 16 
m 



PA; 



12 



= ln(fcA;') - ln(fc„fc;) + 1 In 2, 

61 — a 61 — a 6 

which in turns yields the following final expression for the Renyi entropy. 



(30) 



^^^'^^ = 6 ^ ^ ' 3 1^ ^1(0, J ■ ^ ^ 



Here, the elliptic parameter k is defined in ([7]), fc' = Vl — A;^, the modulus parameter q 
is given by equation (1201) where /(/c) is the complete elliptic integral (E]), and the theta 
functions 9j{z,q) are defined by the series (1221- [25]) . 
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2.2. h<2 



In this case we have 



Am = tanh(m7rro) 



^zme _ I 



where, as usual, 



e = TTTo. 



(32) 



(33) 



The entropy is: 



S 



R 



y In 



m=—oo 
oo 



1-a ^ 



m=l 

oo 



1 + e-2™^ 
1 



+ 



-2me \ 



1 +e- 



-2me 



+ 



1 + e-^"^^ 



1 + e- 

oo 



-2ine 



1-a V 2° 



^ 5: In [1 + e--] - -i^ 5^ In [1 + e-"] + In 2. (34) 



m=l m=l 

Again the second term can be immediately summed using (|T2l) : 



oo oo 



m=l 



m=l 



1 a 



6 1 — a 



In 



Ing 



where, as usual. 



q = e 



-■Kl{k')/I{k) 



(35) 



(36) 



The first term, as in the previous case, admits the similar representation involving the 
elliptic parameter ka = k{q°'), 



y In h + e-2"-^l = 

m=l 

Using (l35l) and (1371) in (1341) . we obtained that 
5« = l^ln^ fc^ \ 1 « 



In 



a In g 



(37) 



ln( 1 +ln2 



61 -a yiQk'^J 6 1-a VISA;' 

1 a 11 fkl\ 1 

n — + n + - n 2, 

61-a Vfc2y^61-a Vfc'/ 3 



(3^ 
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which in turns yields the following final expression for the Renyi entropy in the case 
h<2. 

a) = g — In j + - — In ^^^^jj^^^ j + 5 In 2. (39) 

Here, as before, the elliptic parameter k is defined in ([7]), /c' = vT^^X^, the modulus 
parameter q is given by equation (1201) where I{k) is the complete elliptic integral (JHl), 
and the theta functions Oj{z\q) are defined by the series (l22l - 1251). 

Remark. One can wonder about an apparent tautological character of the 
formulae ( l3Ti) and ( l39l) . Indeed, they seem just to re-express one g-series (Sn^pAjd)) 
in terms of the another {6j{0\q)). The important point however is that the g-series 
representing the theta-constants place the object of interest, i.e. the Renyi entropy, in 
the well - developed realm of classical elliptic functions. In fact, to solve a problem in 
terms of Jacobi theta- function is as good as to solve it in terms of, say, elementary 
exponential function (after all, the exponential function is also an infinite series !). 
The crucial thing is that a lot is known about the properties of the theta-constants and 
this allows a quite comprehensive study of the Renyi entropy both numerically and 
analytically. In the next section we will demonstrate the efficiency of equations ( 1311) 
and (|39|). 

3. Renyi Entropy. The Analysis. 

When studying the analytic properties of the Renyi entropy with respect to the variable 
a, it is convenient to pass from the modulus parameter q to the (more standard) 
modulus parameter r defined by the relations, 

g = e--, r = i^-j^ = iro, > 0. (40) 
I{k) 

The theta functions 6j {z, q) then become the functions, 

e^{z\T):=e^{z,e^'^), J = 1,2,3,4, (41) 
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which are holomorphic for all z and for all r from the upper half plane, 

> 0. (42) 

Using these new notations, the above obtain formulae for the Renyi entropy can be 
rewritten as 

for /i > 2 and 

ei(0|mro) \ 1, , ,,,, 

SR{pA,a) = - In - + - In " . + -ln2, 44 

6 1 — a / 3 1 — a \ 6^3 (0|a;zro) 6^4(0 |a«ro)/ 3 

for h < 2. To proceed with the analysis of these expressions as functions of a, we will 
need some pieces of the general theory of Jacobi functions 6j{z\T) which we collect in 
Appendix A. 

Our first observation is that the domain of analyticity (H2|) and the positiveness of 
the parameter tq indicate that the all three theta-constants, i.e. 92{0\aiTo), 63{0\aiTo), 
and 6'4(0|airo) are analytic in the right half plane of the complex a - plane: 

> 0. (45) 

Simultaneously we notice that inequality ( 1A.23I1 implies that the theta-ratios appearing 
in the right hand sides of (H3|) and (jHj) are never zero. Hence we can claim that the 
Renyi entropy, as a function of a, is analytic in the right half plane (jlS]), with the 
possible pole at a = 1. However, since as a 1 the theta-ratios in fH3l) and (jH]) 
become the square roots of the product kk' and of the ratio k'^/k', respectively (see 
also ( l30l) and (|38ll). the singularity at a = 1 is, in fact, removable and we can write 
that 

SniPA, 1) = Hkk') + 1 In 2 + i A In ^^^(Ol^^-o) ^4(0|mro) 



6 ^ ' ' 3 ' 3da V 9l(0\aiTQ 



-Un{kk') + ^ln2 + lAin(A;X)U=i (46) 
o 3 baa 
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for /i > 2 and 



a=l 



= — In — +-ln2H —In 

6 \k^ 3 6 da \kl 



\a=l 



(47) 



for h < 2. It is an exercise in the theory of elhptic functions to show that the expressions 
on the right hand sides of (H6|) and P7|) are in fact the respective Von Neumann 
entropies calculated in [T9l [33| [7]: 

4 .,..2I{k)I{k') 



SiPA) 
S{pa) 



vr 



> 2 
h<2 



This fact, i.e. the statement that 
lini SR{pA,a) = S{pa) , 



(48) 
(49) 

(50) 



can be of course obtained via much more elementary calculations based on the original 
series representation ( |T^ for Sii{pA,Oi). 

Consider now the two other critical cases: a —>■ and a ^ oo. 



3.1. a 



oo 



The limit of large a is interesting for the single copy entanglement suggested by M. 
Plenio and J. Eisert [35] . In fact, the Renyi entropy contains information about all 
eigenvalues of the density matrix and we can extract the largest eigenvalue [maximum 
probability pm] from the limit a — > oo (S'q(pa) — InpA/)- 

Using the first series from equations (]A.20p - (]A.22p we obtain at once that 
92{0\aiTQ)94{0\aiTQ) 



and 



^i(0|mro) 
0|(O|mro) 



2e-^ (1 + (e-"""«)) 



6*3 (0 1 airo) 6^4 (0 1 aira) 



4e-^ (1 + (e-2"--o)) 



(51) 



(52) 
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as a — > oo, —n/2 < arga < tt/2. Plugging these estimates in and (IHI) and 
recalling that tq = I{k')/ I{k), we arrive at the following description of the Renyi 
entropy in the large a limit. 

a flkk' TT I{k')\ 
S„(p.,a) = — (jln- + -^j (53) 

+0 ( ie-"""° 
Ikk' 71 I(k') ^ ( 1 

TT 77 

a — > OO, < arga < — , 

2 2 

for /i > 2, and 

SH(p..a) = — (glnj^--^j+— ln2 (54) 
6 ""4^ ^ 6 /(A;) ^ V« 

77 77 

a — s> OO, < arga < — , 

2 2 

for h < 2. Alternatively, these estimates can be easily extracted from the original 
series representations, i.e. equations (ITS!) and flM|) . with the help of the identities (fTTj) 
and (fT2|) . In other words, the theta-summation of the series (fT8|) and (IMI) is not really 
needed for the large values of the parameter a. 

Remark. The asymptotic representations fl53l) and fl5^ are only valid for the bulk 
of the XY model, i.e. away from critical lines 7 = or /i = 2. Near the critical points, 
when 7 7^ and /i — 2, or 7 — > and h < 2, the module parameter tq becomes small 
and the estimates (l53l) and (15^ are not valid unless the double scahng condition, 

aro — > 00 (55) 

takes place. 
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3.2. a^O 

This is where the theta-formulae help. Indeed, using the second series from the Jacobi 
identities (1A.20P - (lA.22p . we arrive at the estimates, 
^2 (0 1 airo) ^4 (0 1 airo ) 



2e 1 + O e (56) 



and 



-le-^(l + 0(e-^)) (57) 



^3(0|mro)^4(0|mro) 2 
as a To 0, —7i/2 < arga < 7r/2. These formulae indicate the appearance of a 

singularity of order in the Renyi entropy as a 0. In fact, since we consider the 

limit of a large block of spins, the dimension of the corresponding Hilbert space also 

goes to infinity. This is the reason for which the Renyi entropy has a singularity at 

a = 0. 

Substituting ( !56l) and ( l57l) into ( l43l) and f j44l) . respectively, we obtain the following 

description of the Renyi entropy in the small a limit. 

„ , - 1 vr I(k) a 1 ^ kk' ^ / 2_\ , ^, 

+ '^^^ ' 0(a), (59) 



a 12 1{k') 



TT TT 

a^O, --<arga<-. 



for h > 2, and 



1 + a TT /(/c) . , s 

TT TT 

--<arga<-, 

for h<2. 

Similar to the case of the Von Neumann entropy dealt with in [22], equations (135]) 
and fl60l) can be also used for the evaluation of the small tq = I{k')/I{k) limit of the 
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Renyi entropy with the fixed a > 0. This hmit (cf. [33]) appears either in the case of 
the critical magnetic field, i.e. 7 7^ and /i — > 2, or when approaching the XX model, 
i.e. 7 — i> and h < 2. We shall now consider these limits. 

3. 3. Critical magnetic field: 7 7^ and h 2 

This is included in Case la and Case 2 which means that, 

k = l- ^\h-2\+0(\h-2\^), k' = -\h - 2\^ (1 + 0(\h - 2\)) , (62) 
27"^ 7 

and, in turn, 

= -ln|2-/i| +21n47 + 0(|/i-2|ln^|/i-2|), (63) 

I{k ) 

h^2, 7^0. 

This means that in this limit tq ^ and we can use fl58|l to arrive at the following 
estimates for the Renyi entropy in the case of the critical magnetic field, 

Sr{pa, a) = ^ (-^ In |2 - /i| + 1 In 47) (64) 
a V 12 6 / 

+0{\h-2\ln^\h-2\). 

We notice that the singularity of the Renyi entropy is logarithmic like for the Von 
Neumann entropy, but coefficient in front of the logarithm is different and a-dependent. 

3.4- An approach to XX model: 7 — > and h < 2 
This is included in Case lb which means that, 

= 1-^ + 0(7^), fc' = -^(l + 0(f)), (65) 



and, in turn, 



7r4S: = -21n7 + ln(4-/i2) + 21n2 + 0(7ln2 7), (66) 
I{k') 



7-^0, h<2^/l--f^. 
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Again, since tq — > 0, we can substitute these into (1^ and arrive at the following 
estimates for the Renyi intropy in the case of the XX model limit 

Sr{pa, a) = ^ (-1 In 7 + ^ ln(4 - + i In 2) (67) 
a V 6 12 6 / 

+0(7 In' 7). 



We note that if a = 1 then equations fl64j) and fl67j) transforms to the respective 
formulae for the Neumann entropy obtained earlier in [33]. 

3. 5. The factorizing field 



We already showed in the introduction that for h = hf{'~)) = 2a/1 — 7' the ground 
state can be written as 

\GS) = \GSi) + \GS2) , (68) 

where \GSi^2) are the product states given in and clearly have no en- 
tropy/entanglement by themselves. 

We can calculate the Renyi entropy of the ground state at the factorizing field by 
considering the limit k ^ of ( 14^ . Remembering that, using ( ]A.20tlAl22!) in this limit 



a/2 



^2(0|mr ~ 0) ~ 2 , (69) 

03(O|mr ~ 0) ~ 04(O|mr ~ 0) ~ 1 , (70) 

it is easy to show that 

SniPA, a) = \n2 (71) 
regardless the value of a. This result is not surprising and was to be expected in light 



of (1681) . In fact, the limiting density matrix of the block of spins at the factorizing field 
is (1/2) X I2, where I2 is the 2x2 Identical matrix. 

Please note the importance of the order of limits around the factorizing field. In 
fact, the expression in flTTl) is independent of a and therefore regular in the limit a —>■ 0, 
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while off the factorizing field line the entropy diverges like in (1611) for a — >■ 0. As one 
approaches the factorizing field, k ^ and therefore Tq ^ oo in such a way that aro 
stays constant. 

4. Renyi Entropy and the Modular Functions. 

The square of the elliptic parameter k, considered as a function of the modulus r, is 
usually dented as A(r), and it is called the elliptic lambda function or A - modular 
function. We note that (cf. ( I2T1) ) 

A(r) = = fc^(e-^), -sr>0, (72) 

and that 

l-AW^||t|.*'^e-). (73) 

The function A(r), sometimes also denoted as k^(t), plays a central role in the theory of 
modular functions and modular forms, and a vast literature is devoted to this function 
- see the classical monograph [36]; see also [31], [37], [38], [39] and Section 3.4 of 
Chapter 7 in [40j. The function A(r) possess several remarkable analytic and arithmetic 
properties, some of which are listed in Appendix B. 

In terms of the A - modular function, the formulae for Renyi read as follows. 

SniPA^a) = ^ ln(fc k') - ^ ln(A(mro)(l - A(mro))) +]-\n2, (74) 

61 — a 12 1 — aV /3 

ioT h > 2 and 

r,/ \ I a ^ f k'\ 1 1, A^(mro) 1, ^ 

SR{pA,a) = -- In +—- In- ^-^ + -ln2, 75 

61 — a \k J 12 1 — a 1 — X{mTQ) 3 

for /i < 2. These relations allow to apply to the study of the Renyi entropy the 
apparatus of the theory of modular functions. We are going to address this question 
specifically in the next publications. Here, we will only present the two most direct 
applications of the modular functions theory related to the symmetry properties of the 
A-function indicated in (m - flRTl) ). 
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4.1. Modular transformations 
Put 

/(r):=A(r)(l-A(r)), and 9{r) = (76) 

and re-write the formulae for the Renyi entropy one more time: 

SRipA.a) = ^ ln(fc k') - ^ ln/(mro) + ^ln2 (77) 

b 1 — a 1/1 — a 6 

for h > 2, and 

,la/A;'\ll,.,l, 

SRipA^a) = ^ -. 7^ + 777 1 In^(mro) + -ln2 (78) 

b 1 — a \/c^/ 12 1 — a o 



for h < 2. The symmetries flB.2p and fIB.Sp imply the following symmetry relations for 



/(r) and g{T) with respect to the action of the modular group, 



/ ^--J = fir) (80) 
(7(r + l) = (7(r), (81) 

It follows then, that the function /(r) is automorphic with respect to the subgroup of 
the modular group generated by the transformations, 

T^T + 2 and r ^ --, (83) 

r 

while the function g{T) is automorphic with respect to the subgroup of the modular 
group generated by the transformations, 

T ^ T+1 and T — - — . (84) 

2r + 1 ^ ^ 

Of course, the both functions inherit from the lambda-function the automorphicity 
with respect to subgroup (1B.6I) (which is a common subgroup of the subgroups (183|) 
and flM|) ). Therefore, we arrive at the following conclusion. 
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Proposition. Up to the trivial addition terms and multiplicative factors, and 
after a simple re-scaling, the Renyi entropy, as a function of a, is an automorphic 
function with respect to subgroup (E^i of the modular group, in the case h > 2, and it 
is automorphic with respect to subgroup [84\ ) of the modular group, in the case h < 2; 
in both cases the entropy is automorphic with respect to subgroup W.6\) . 



The indicated symmetry properties of the Renyi entropy yield, in particular, the 
following explicit relation between the values of the entropy at points a and l/aT^. 

Sr[pa,— \= / ri - « Sfi PA, « + ^ In -r. 85 

\ oTq / arg — 1 6 oTq — 1 4 



for h > 2 and 



"^^ln/(mro). 



12 - 1 



for /i < 2. We bring the attention of the reader to the appearance in the case h <2 oi 
an extra term involving the modular function /(r) . 

4.2. a = 2" 

For the indicated values of the parameter a one can apply Landen's transformation 
(1B.7P and reduce the function A(airo) to the function 



Hence, for these values of a the Renyi entropy becomes an elementary function of the 
initial physical parameters h and 7. Let us demonstrate this in the case a = 2. 
From (]B.7p it follows that 



\{2lTa) 



1 + k' 
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Therefore, 



/(2«ro) 



4k'{l - k'f 
{l + k'Y 



and g{2iT()) 



Ak'{l + k'Y' 




Using these, we can find the Renyi entropy for a = 2 from fl77|) and (178|) : 




(87) 



for /i > 2 and 



Snia = 2) = -- In 



) 



+ -ln2 
2 



8. 



;8 



k^ (1 + A;') 



for /i < 2. Repeating Landen's transformation again and again, we can iteratively 
construct a ladder of "elementary" entropies for increasing values of a = 2". 

5. Summary and Conclusions 

We analyzed the entanglement of the ground state of the infinite one-dimensional XY 
spin chain by calculating the Renyi entropy Sci{pa) of a large block A of neighboring 
spins. The Renyi entropy has been proposed as a meaningful measure of the quantum 
entanglement of a system and it is a natural generalization of the Von Neumann entropy. 
In fact, for a = 1 the quantities are equal. Moreover, knowledge of the Renyi entropy 
for all a's allow for the reconstruction of the density matrix and an easier identification 
of the sources of entanglement in the mixed quantum state. 

We arrived at an analytic expression of the entropy in the bulk of the two- 
dimensional phase diagram of the model, in terms of an elliptic parameter and elliptic 
theta functions. These expressions allowed us to study the behavior of the Renyi 
entropy for the different values of a and of the parameters of the model. We found 
the limiting behavior of the entropy for a ^ oo, which is essentially the single copy 
entanglement introduced in [35j. In that work, it was shown that this quantity scales 
like 1/6 In L for the isotropic XX model. This is consistent with our findings - setting 
~ 1, fc' ~ in ( !54l) - and we generalize it to the rest of the phase diagram. 
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In the hmit a — >■ we showed that the entropy diverges hke a~^. A very interesting 
behavior occurs at the factorizing field /i/(7) = 2a/1 — 7^. On this line the ground 
state can be written as a sum of two product states. This means that the reduced 
density matrix remains proportional to the two-dimensional identity matrix and we 
showed that the Renyi entropy is 5*^ = In 2, independent from a. So, even for a ^ 
the Renyi entropy stays finite at the factorizing field, while it diverges as one moves 
away from this line. 

The bulk of the XY model is gapped and the entropy of a large block is known to 
saturate to a finite value, which we calculated. As one approaches the critical lines, the 
entropy diverges logarithmically in the gap size. We calculated exactly the prefactor 
of this logarithmic divergence as a function of a for the two universality classes of the 
critical lines and found agreement with the Von Neumann result at a = 1, as to be 
expected. 

Finally, using the properties of the theta functions, we showed that the limiting 
Renyi entropy is a modular function of a. The properties of the entropy under modular 
transformations seem very interesting and will be the subject of a subsequent paper. 
In a previous work [23] we showed that the curves of constant entropy are ellipses and 
hyperbolae and that they all meet at the point {h,'y) = (2,0), which is a point of 
high singularity for the entropy. This is valid also for the Renyi entropy and seems 
to be connected with the aforementioned modular properties of the entropy. We will 
investigate this relationship in a future work. 
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Appendix A. Theta Functions 

In this appendix the necessary facts of the theory of Jacobi theta-functions are 
presented. For more detail, we refer the reader to any standard text book on elliptic 
functions, e.g.[3l]. 

Among the four theta-functions, only one is functionally independent, and usually 
it is taken to be the function 63{z\t). The rest of the theta-functions are related to 



6*3 (2; I r) via the simple equations, 

e,iz\T) = -ze"f + 1 + ^ TTT , (A.l) 

d2{z\T) = e^+'^e, + ^ vrr , (A.2) 

9,{z\T) = 9s(^z + ^n\Ty (A.3) 

The principal characteristic properties of the theta-functions are their quasi - 
periodicity properties with respect to the shifts, z z + n and z z + vrr: 

9,{z + 7i\t) = -9i{z\t), (A.4) 

9i{z + 7rr|r) = -e-^'^-^"9i{z\T), (A.5) 

92{z + 7r\r) = -92iz\T), (A.6) 

92iz + 7rr|r) = e-'''^-^''92iz\r), (A.7) 

^3(^ + 7r|r) =^3(^|r), (A.8) 

9siz + nT\T)=e--'^-^''9s{z\T), (A.9) 

94{z + n\T) =94{z\t), (A.IO) 
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e^iz + 7tt\t) = -e-"'"-"*"04(2|r). (A.ll) 

The complementary set of the properties is the set of the following symmetry relations 
with respect to the transformations, r ^ r + 1 and r —t~^ (that is, with respect 
to the action of the modular group): 

(A.12) 

^ '^■e'^ei{z\r), (A.13) 





+ 1) = 










z\t 


+ 1) = 











t \ t 

eTe2{z\T), (A.14) 
02[^^\-^)=^e'^9,iz\r), (A.15) 

es{z\T + i) = e^{z\T), (A.16) 

ye-^3(z|r), (A.17) 

(A.18) 





-7) 






z\t 


+ 1) = 











,2 



-e-^2(2|r), (A.19) 



where the branch of the square root is fixed by the condition. 



= 1, if T = i. 



An immediate important corollary of these relations is the possibility of the following 
alternative series representations (the Jacobi identities) for the theta-functions 
participating in the formulae (j^3l) and (jUj) for the Renyi entropy. 



9,{0\t) = 2 e^^^("+^)' = JUl + 2 ^(-l)"e-"^ ) , (A.20) 

n=0 \ n=l 

oo r~r / oo 

^^3(0|r) = 1 + 2 5^ e— ^ = V M ^ + 2 S """^ ) ' (^-^^^ 

n=l \ n=l 

oo r~r oo 

,2 „ / z 



, T 

n=l n=0 



^4(0|r) = l + 2^(-l)"e"'™' = 2W-^e-^("+2) . (A.22 
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The first series in each of these formulae allow an efficient evaluation of the 
corresponding theta-constant for large Qt, while the second series provides a tool for 
analysis of the theta-constant in the limit of small |r|. In Section 4 we use these 
identities for investigating the singularity of the Renyi entropy at a = 0. 

The last general fact of the theory of elliptic theta-function we will need, is the 
description of their zeros, as the functions of the first argument. In view of the relations 
( 1A.1|) - (]A.3|) it is sufficient to describe the zeros of 93{z\t). They are: 



z = Znm = - 7r + -7rr + n7r + mvrr, n, m (H5i) G Z. 



This information about the zeros of 63{z\t), taking in conjunction with the relations 
( ]A.2I) and (1A.3I) . implies, in particular, that 

^^2(O|r)03(O|r)^^4(O|r) ^ 0, Vr, > 0. (A.23) 
Appendix B. Elliptic Lambda Function 

The properties of the A-function presented below form an important but very far from 
being exhausted set of the extremely exciting properties and connections which this 
function enjoys. For more on the lambda and related functions we refer the reader, in 
addition to the monographs already mentioned, to the websites [41j and [12] and to 
the references and links indicated there. 

(i) Let Q denote the "triangle" on the Lobachevsky upper half r - plane with the 
vertices at the points 0, 1 and oo and with the zero angle at each vertex (the edges 
are: 3?r = 0, 3?r = 1, |r — 1/2| = 1/2). Then, the function w = A(r) performs 
the conformal mapping of the triangle Q onto the upper-half plane '^w > 0, and 
it sends the vertices 0, 1, and oo to the points 1, oo, and 0, respectively. It also 
should be noticed that the real line is the natural boundary for A(r) - the function 
can not be analytically continued beyond it. 
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(ii) The direct corollary of the conformal property just stated is the following analytic 
fact. Let {/, z} denotes the Schwartz derivative, 



f'U) 2 V /'(-- 

Then, the lambda-function A(r) satisfies the following differential equation, 

= "2 - 2 (A^ + A(A^ ^^-'^ 
(iii) The function A(r) satisfies the following symmetry relations with respect to the 

actions of the generators of the modular group (cf. (1A.14P - (1A.19P ). 

A {--\ = 1 - A(r). (B.3) 



These symmetries in turn imply the equations. 



A(T + 2) = A(r), (B.4) 



which show that the function A(r) is automorphic function with respect to the 
subgroup of the modular group generated by the transformations, 

r^r + 2, and — - — . (B.6) 

2r + l ^ ' 

(iv) In addition to the symmetries with respect to the modular group, the function 
A(r) satisfies the so-called second order transformation, also called Landen's 
transformation, which describes the action on A(r) of the doubling map, r 2r: 



Here, the branches of the square roots are chosen according to the equations (cf. 
1 and (CI), 



V^=^-Me ), Vl-A(r) = ^.fc(e 
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(v) By means of the algebraic equation, 

T(r) - ^ (l-Mr) + AW 

^^^^ - 27 A2(r)(l-A(r))2 ' ^^"^^ 
the elhptic lambda - function defines even more fundamental object of the theory 
of modular forms - Klein's absolute invariant J{t). The function J(r) is a modular 
function, i.e. it is automorphic with respect to the modular group itself, 

J(r + l) = J(r), j(^-i^=J(r); (B.9) 

moreover, any other modular function is algebraically expressible in terms of the 
invariant J{t). The function J(r) admits also an alternative representation in 
terms of the Ramanujan-Eisenstein series Ej: 

We remind that 

oo oo 

E,{t) = 1 + 240 Y,(^3{k)q'^ E,{t) = 1 - 504 ^ a5(A;)g2^ q = e'^\ 

k=l k=l 

where (Tk{n) is a divisor function, i.e. 
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